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Abstract. We present radial velocity and new interferometric measurements of the double star Atlas, which permit, with the 
addition of published interferometric data, to precisely derive the orbital parameters of the binary system and the masses of the 
components. The derived semi-major axis, compared with its measured angular size, allows to determine a distance to Atlas of 
132 + 4 pc in a purely geometrical way. Under the assumption that the location of Atlas is representative of the average distance 
of the cluster, we confirm the distance value generally obtained through main sequence fitting, in contradiction with the early 
Hipparcos result (1 18.3+3.5 pc). 



Key words. Stars: individual: Atlas 
tions: individual: Pleiades 



Stars: distances - Stars: fundamental parameters - Galaxy: Open clusters and associa- 



1. Introduction 

The distance of the Pleiades has been under debate ever since 
the publication of Hipparcos parallaxes (ESA 1997). Hipparcos 
data allowed to precisely calibrate and compare the posi- 
tions of open cluster main sequences with those obtained us- 
ing the usual main sequence fitting technique, and resulted 
in a mismatch in the case of at least two clusters (Pleiades 
and Coma Ber), the mos t severe being that of the Pleiades 
( Pinsonn eault et al1ll998l). Hipparco s locates the Pleiades at 
118.3 ± 3.5 pc Jvan Leeuwenl fl999). abo ut 10% nearer tha n 
usually quoted results (132+4 pc, see e.g. lMevnet et al.ll9 93 ), 
that is about 14 pc nearer, or ~ 0.3 mag brighter in distance 
modulus, or ~ 1 mas larger in parallax. 

Various explanations have be en proposed to solve th e prob- 
lem: large helium abunda nce jMermilliod et alj|l997h . metal 
deficiency (Grenon 2000), problem in Hipparcos parallaxes 
(Makarov 2002). None of them brought scientific consensus. 
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Th e elements in the balance are not minor ones, as pointed 
out by IPaczvnskl J2003I) . On one side, stellar evolution the- 
ory and photometric calibrations are at stake, or, on the other 
side, problems appear when averaging parallaxes in order to 
reach sub milliarcsecond accuracy from Hipparcos data. The 
problem must be solved, especially in the perspective of the 
future Gaia mission, which will acquire astrometric data with 
the same principle as Hipparcos. Thus, identifying the source 
of the problem is worth the effort. It is perhaps relevant to note 
here that the Pleiades lie near the ecliptic, where the number of 
great circles abscissae used for astrometry is especially small. 

Fortunately, there exists several other ways to determine 
distances with no or little model dependencies. 

This paper presents radial velocity (hereafter RV) measure- 
ments of the interferometric binary Atlas, which permit to de- 
rive the distance of Atlas from purely geo metric arguments . 
This star was first discovered to be a binary by B artholdJ d 1 975h 
from occultation by the Moon. More recent interferometric 
measurements bv lPan et al J (12004) yielded well constrained or- 
bital parameters which lead to a distance of 135 + 2 pc, but 
only when combined with a mass-luminosity relation, making 
the result model dependant. On the contrary, our RV data makes 
our work completely model free. 

Another star allowed an independent measurement of the 
dist ance to the P leiades. Using Hipparcos epoch photome- 
try, iTorresI J2003I) discovered that the spectroscopic binary 
HD 23642 (HIP 17704) is also eclipsing. iMunari et al J ( 12004 
obtained photometric and additional RV observations and de- 
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Table 1. New astrometric measurements taken at the Markm 
(BJD < 2450000) andNPOI (BJD > 2450000) interferometers. 
BJD means Bary centric Julian Day. Column 2 gives the number 
of visibilities, while the last three columns give the axes and 
orientation of the error ellipse. 
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rived a distance of 132 + 2 pc, perfectly compatible with the 
result of the main sequence fitting technique, even though the 
error bar looks rather optimistic. This method is almost model- 
free, but not entirely since the effective temperatures of the 
components were determined through comparison with syn- 
thetic spectra. 

In addition, Soderblom et al. recently announced (AAS 
meeting) a distance of 135 + 3 pc on the basis of new paral- 
lax measurements made with the HST on three stars. 

2. Observations and data reduction 

2.1. Additional interferometric observations: 

Atlas was observed by CH between 1989 and 1992 with the 
Mark in interferometer, and in 1997 and 2000 with the Navy 
Prototype Optical Interferometer (NPOI). The journal of these 
astrometric observations is given in Tabled The reduction and 
calibrati on of t he dat a followed the procedures described by 
iHummel et alJ Jl998l) . A lthough these da ta are less numerous 
than those published by IPan et al] J2004I) on the basis of ob- 
servations with the Mark in stellar interferometer and with the 
Palomar Testbed Interferometer (PTI), those collected with the 
NPOI have the advantage of being free from any ambiguity 
on the orientation of the orbit, because NPOI measures clo- 
sure phases. The precision is comparable with that reported by 
Pan et al . (2004), so we could merge both samples together and 
obtain an improved estimate of orbital parameters (see Table 0] 
below), using the OYSTER code. The observations and fitted 
orbit are shown in Fig.^ 

2.2. Spectroscopic observations: 

Echelle spectra of Atlas were obtained using above all the 
ELODIE spectrograph attached to the 1.93m telescope at 
Observatoire de Haute-Provence, France, but also with the 
CORALIE spectrograph attached to the Swiss Euler 1.2m tele- 
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Fig. 1. Astrometric orbit with published data from IPan et alJ 
(2004) (dotted error ellipses) and new data from C. Hummel 
(full lines; see Table[0. The ellipses indicate the errors on each 
axis, while the full straight line denotes the periastron. 



scope at ESO La Silla, Chile, and with the FEROS spectrograph 
attached to the 2.2m telescope at ESO La Silla. Almost all ob- 
servations were gathered during the 2003-2004 season, though 
a few archive spectra taken with ELODIE date as far back as 
1995. All spectra were taken without simultaneous ThAr cali- 
bration, either in star mode only, or in star+sky mode. Until the 
end of February 2004, observations were made within time nor- 
mally allocated to other programmes, thanks to the good will 
of some observers. From the 27th of February to the 29th of 
March 2004, spectra were taken at OHP within officially allo- 
cated time in service mode. The exposure times were typically 
10 - 15 mn, resulting in a S/N ratio varying from ~ 120 in the 
worst atmospheric conditions to more than 400 in the best ones. 

The ELODIE and CORALIE spectra were reduced and nor- 
mali zed to con tinuum using IRAF, with the method described 
bv lErspamer & North (2002). The FEROS spectra were re- 
duced using the standard reduction pipeline provided by ESO 
in the MIDAS context. Then, a more careful continuum nor- 
malization was performed in the immediate vicinity of the few 
lines which appeared useful for RV determination. 

The journal of the observations is given in Table |2] Since 
spectra with S /N <~ 200 proved to be of little help because of 
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Table 2. Journal of the specroscopic observations. Epochs 
are given in Barycentric Julian Day. Key to the observers: 
yd=Y. Debernardi, cv=C. Vuissoz, fg=F. Galland, mg=M. 
Groenewegen, rm=R. Monier, pn= P. North, rds=R. da Silva, 
bp=B. Pernier, pg=P. Girard 
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Table 3. Lines used in this work for radial velocity estimations. 
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the shallowness of the lines (< 0.1), only the best exposed ones 
are listed here. 

3. Analysis of spectra and orbit determination 

Since the lines of the components are wide and alw ays blended, 
we have used the KOREL code by Hadrava ( 1995) for spectral 
disentangling and RV determination. Even with this tool, the 
task was a delicate one because of the very small number of 
suitable lines and because of their shallowness. Unfortunately, 
the Balmer lines could not be used because of their larger width 
compared with the relatively small RV amplitude. The lines 
we used are listed in Table [3] The cleanest is that of Mg n. 
The two Si n lines are appropriate as well, but suffer from con- 
tamination by telluric lines which, although not very deep, do 
hamper convergence if they are not corrected for. In order to 
suppress these telluric lines, we used the high resolution spec- 
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Fig. 2. Disentangled lines of Mg n centered on the systemic ve- 
locity. Full line: primary; broken line: secondary. Dotted line: 
spectrum observed in 1995 near quadrature, at phase ~ 0.1. 

trum of telluric lines provided by Hinkl e et alJ (2000) together 
with their visible atlas of the Arcturus spectrum. This spec- 
trum was convolved with a gaussian instrumental profile corre- 
sponding to the resolving power of ELODIE R = 42000; since 
CORALIE and FEROS have only slightly better resolving pow- 
ers (R = 50000 and 48000 respectively), the same convolved 
telluric spectrum proved useful for them as well. Since airmass 
and humidity varies from one spectrum to the next, the correc- 
tion had to be adjusted visuall y for each spectru m by scaling 
the telluric line intensities (see iNorfh e t al. 1994, for details), 
mainly on lines lying on the continuum. This procedure proved 
quite satisfactory, thanks to the fact that the relative depths of 
the telluric lines remain the same in a given spectrum even 
though their absolute depths change from one spectrum to the 
other. 

Fortunately, most orbital parameters are known from inter- 
ferometry, especially To, e, i and co, so these quantities could 
be fixed while running KOREL, leaving only the RV amplitude 
K\ and the mass ratio q as adjustable parameters. 

By limiting the data set to spectra with S /N > 200, and 
with starting values of K\ between 10 and 60 km s 1 and of q 
between 0.3 and 1.0, we found good convergence towards the 
values given in Table 0] The disentangled spectra are shown 
in Figure |2 for the Mgn line. One clearly sees that the pri- 
mary rotates much faster than the secondary, the line exhibit- 
ing a typical rotational profile. The vsinz values are respec- 
tively about 240 and 60 kms -1 . The RV curves computed by 
the KOREL code are shown in Figure|5J the points are obtained 
by cross-correlating the individual spectra with the disentan- 
gled spectra. The error estimate on the fitted orbital parameters 
is not straightforward, since the disentangling method used in 
the KOREL code does not give it in a nat ural way. Therefore, 
we ha d recourse to the jackknife method (Efron & Tibs hiranil 
1993), which consists in removing each spectrum in turn from 
the whole set of data, and run KOREL with the remaining 23 
spectra. One obtains slightly different values of the fitted pa- 
rameters, and their variance multiplie d by the nu mber of spec- 
tra gives the final error estimate (e.g. Iliiic 2003). 



Date BJD Instr. expos. S/N obs. g 

-2400000 time(s) | 0.95 
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Fig. 3. Radial velocity curves given by the KOREL code. All 
orbital parameters were fixed at their astrometric values, except 
K\ and q. There are three radial velocity estimates per epoch, 
one per line, even though the RV curves were constrained si- 
multaneously by all three lines. Note that the systemic velocity 
is arbitrarily set to zero because the KOREL code cannot deter- 
mine it. 

Table 4. Astrometric and spectroscopic orbital parameters ad- 
justed with the OYSTER and KOREL codes respectively. A 
simultaneous fit of astrometric and RV data yields essentially 
the same results. The errors on the spectroscopic parameters 
are those obtained through the jackknife method alone, but the 
error on d includes the small uncertainty on a". 
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Parameter 
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Period [d] 


290.984 ± 0.079 


K, (kms^ 1 ) 


26.55 ± 1.41 


To [BJD] 


2450583.0 ± 1.9 


K 2 (kms -1 ) 


36.89 ± 0.22 


e 
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154.0 ± 0.7 


M 2 (M e ) 


3.42 ± 0.25 


a" (mas) 


13.08 ±0.12 


d(pc) 


132 ±4 



Since KOREL gives individua l RV values, we mixed the 
latter with the astrometric data of P an et alJ J2004I) and ours 
to fit simultaneously all orbital parameters. This is in princi- 
ple a more elegant method than to fit separately the spectro- 
scopic data and the astrometric ones. However, in our case 
one has to bear in mind that the individual RV values are not 
mutually independent stricto sensu, because they result from 
a cross-correlation between each individual spectrum and the 
components' disentangled spectra. Since the latter are built on 
the same data, there is some circularity which tends to reduce 
the O - C residuals. In addition, the spectroscopic orbit is not 
independent from tbe astrometric one, since the latter was used 
to constrain most orbital parameters. Nevertheless, such a con- 
sistency check is reassuring in that it provides the same result 
as when combining the separate astrometric and spectroscopic 
solutions. 



4. Distance of Atlas and conclusions 

The orbital parameters obtained above yield a final distance 
of d = 132 + 4 pc which confirms the "traditional" distance 
obtained through t he main s equence fitting technique. Thus the 
result obtained bv lPan et"al] J2004I) is fully confirmed, and there 
is no need for any revision of the stellar evolution models, or 
assumed helium content of the Pleiades members, or to postu- 
late a non-solar metallicity. 

Our result has the advantage of relying exclusively on em- 
pirical data, without the least recourse to any model. The only 
assumption we made is that of no intrinsic spectral variability 
of the components of Atlas, an assumption which is well jus- 
tified a posteriori for the Mg n line at least, but which will be 
examined further in a more complete forthcoming paper. 

It is worth mentioning that our result is completely coher- 
ent with the distance of the recen tly discovered eclipsing binary 
HD 23642 dMunari et al.ll2004T) . within the errors and within 
the expected depth of the Pleiades cluster. Thus, there are now 
two objects for which a more or less direct distance estimate 
puts the Pleiades at about 130 - 135 pc, which is clearl y dif- 
feren t from the Hipparcos estimate published by Ivan Leeuwenl 
dl999h . It appears increasingly urgent to understand what went 
wrong with that est imate, a difficult task which is in progress 
(Ivan Leeuwenl20 04 ) . 
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